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Abstract. 

We study the ground state phase diagram of a two-leg spin ladder with anisotropic 
ferromagnetic leg couplings under the influence of a symmetry breaking transverse 
magnetic field by the exact diagonalization technique. In the case of antifcrromagnetic 
coupling between legs we identified two phase transitions in the plane of magnetic 
field vs interchain coupling strength. The first one corresponds to the transition from 
the gapped rung-singlet phase to the gapped stripe-ferromagnetic phase. The second 
one represents the transition from the gapped stripe-ferromagnetic phase into the fully 
polarized ferromagnetic phase. 



PACS numbers: 75.10.Jm Quantized spin models, 75.10.Pq Spin chain models 
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1. Introduction 

There has been recently considerable interest in study of the magnetic field-induced 
effects in low- dimensional quantum spin systems. In particular, critical properties of 
the spin S — 1/2 isotropic antiferromagnetic two-leg ladders in a magnetic field have 
been a field of intense studies. This seems pertinent in the face of great progress made 
within the last years in fabrication of such ladder compounds pp. Moreover, since the 
antiferromagnetic two-legs ladder systems have a gap in the spin excitation spectrum, 
they reveal extremely rich quantum behavior in the presence of magnetic field (See 
for recent review [2]). Such quantum phase transitions in spin systems with gapped 
excitation spectrum were indeed studied experimentally [3-8]. On the theoretical side 
these transitions were intensively investigated using different analytical and numerical 
techniques [9-23]. 

Ladder systems with ferromagnetic legs are much less studied. Partly, the reason 
is that spin-ladders with ferromagnetic legs are not still fabricated. However, from the 
theoretical point of view these systems are extremely interesting, since they open a 
new wide polygon for the study of complicated quantum behavior, unsuspected in more 
conventional spin systems [9,24-30]. The variety of open possibilities is clearly seen from 
the weak-coupling phase diagram of a two-leg ferromagnetic ladder which is a function 
of the intraleg exchange anisotropy (A) and the isotropic interleg coupling (JjJ [28], it 
is presented in FigJTJ In addition to the fully gapped rung-singlet phase (RS) (common 
in the case of antiferromagnetic ladder [HIED]), the ground state phase diagram contains 
the gapless Spin-Luttinger-Liquid phase with the easy-plane anisotropy (SLL) and the 
stripe-ferromagnetic (SFM) phase which are realized only in the case of ferromagnetic 
legs (A > 0)0 EE]. 

The very rich and interesting new phenomena arise in the case of competing 
ferromagnetic intraleg and antiferromagnetic interleg couplings, in particular in the 
presence of symmetry breaking magnetic field [28l [29] . The effect of the uniform 
magnetic field, applied parallel to Z (quantization) axes, on the properties of the two-leg 
ladder systems was first studied (in the case of equivalent spin 5" = 1 Heisenberg chain 
model) by H. Schulz [9]. In the case of gapped rung-singlet phase the magnetization 
appears only at finite critical value of the magnetic field equal to the spin gap [9j [12] . 
This behavior is generic for the spin gapped U(l) symmetric systems in a magnetic 
field which leaves the in-plane rotational invariance unchanged [30] and belongs to the 
universality class of the commensurate-incommensurate (C-IC) transitions [3"T] . 

The effect of an uniform transverse field in the case of U(l) symmetric phase is 
highly nontrivial. In the case of classical spin chains this effect has been studied already 
a decade ago [32]. However, in the case of quantum antiferromagnetic XX Z chain this 
problem is still the subject of intense recent studies [32-41]. 

In this paper, we address the similar problem and study the effect of uniform 
transverse magnetic field on the ground state phase diagram of a two-leg ladder with 
anisotropic, ferromagnetically interacting legs coupled by antiferromagnetic interleg 



Spin Ladder with Anisotropic Ferromagnetic Legs in a Transverse Magnetic Field 3 




Figure 1. Schematic picture of the ground state phase diagram of the two- leg ladder 
as a function of the intraleg exchange anisotropy (A) and isotropic interleg coupling 
( Jj_) in the considered in this paper sector of the phase diagram corresponding to the 
competing ferromagnetic intraleg (A > 0) and interleg antiferromagnetic (Jj_ > 0) 
exchange. The Spin-Luttinger-Liquid phase is denoted by SLL, while the Stripe- 
Ferromagnetic phase - by SFM. (From Ref [28].) 

exchange. The Hamiltonian of the model under consideration is given by 

H = < } + < } + H ± C 1 ) 
where the Hamiltonian for leg "j" is 

N 

H-leg = ~ [ ( S i,n S i,n+l S j,n S j,n+l + ^ S j,n S j,n+l) + h^^.n] (2) 

n=l 

and the interleg coupling is 

N 

ffj. = '±EMX* + <A + <^) • ( 3 ) 

Here sj** are spin s = 1/2 operators on the n-th rung, and the index j = 1,2 denotes 
the ladder legs. The intraleg coupling constant is ferromagnetic, J > 0, and therefore 
the limiting case of isotropic ferromagnetic legs corresponds to A = 1 and the external 
magnetic field is proportional to h = Jh ext . 

The model ([1]) has been studied recently using the continuum-limit bosonization 
approach [28j [29]. The main attention in this studies was focused on the investigation 
of new field induced effects in the sector of phase diagram corresponding to the gapless 
easy-plane phase, present in the case of weak antiferromagnetic intraleg exchange 
< Ji < J P [28] . It has been shown, that in the presence of an infinitesimally small 
transverse magnetic field which breaks the in-plane rotational symmetry the gapless 
phase is unstable towards the gapped stripe-ferromagnetic (SFM) phase [29]. The SFM 
phase is characterized by the uniform magnetization in the direction of applied field 
and opposite magnetization of legs in the in-plane direction perpendicular to the field. 
It has been also shown, that when the magnetic field exceeds some critical value, the 
interleg antiferromagnetic order disappear and the system passes into the fully polarized 
ferromagnetic phase [29] . 
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In the opposite case of strong rung exchange Jj_ 3> J an analytical description 
of the magnetic properties of the system in the presence of transverse magnetic field 
is available only in the SU(2) symmetric case at A = 1. In this limiting case it has 
been shown that the isotropic ladder in a magnetic field shows two second order phase 
transitions: at h ext = h^f* from a spin gapped rung-singlet (RS) phase to a gapless 
Spin-Luttinger-Liquid (SLL) phase and at h ext = h^f* a transition from a SLL phase 
into the fully polarized ferromagnetic (FM) phase (see Fig 12 a) [29J . 

In the case of finite intraleg exchange anisotropy (A ^ 1), when the effect of 
symmetry breaking transverse field is important, an analytical solution is not available. 
However, based on the qualitative estimations and symmetry analysis it has been shown 
that the two-step transition from the RS into the FM phase is present also in the case 
of ladder with anisotropic legs [29]. In the same work it has been also predicted that 
for A ^ 1 the RS phase is separated from the FM phase by the gapped SFM phase (see 



In this paper we continue our studies of the anisotropic spin ladder in the transverse 
magnetic field. In particular we apply the modified Lanczos method to diagonalize 
numerically finite (N = 12, 16, 20, 24) ladder systems. Using the exact diagonalization 
results, we calculate the spin gap, magnetization, various spin-structure factors and the 
rung dimerization order parameter as a function of applied transverse field for various 
values of the anisotropy parameter A / 1. Based on the exact diagonalization results 
we obtain the magnetic phase diagram of the ferromagnetic ladders with anisotropic 
legs showing two phase transitions in the plane of magnetic field vs interchain coupling 
in agreement with the predictions made in Ref. [27]. 

The outline of the paper is as follows: In section II we briefly discuss the model in 
the strong rung coupling limit and derive the effective spin chain Hamiltonian to outline 
the symmetry aspects of the problem. In section III we present the results of exact 
diagonalization calculations using the modified Lanczos method. Finally we conclude 
and summarize our results in section IV. 

2. Large rung coupling results 

In this section we briefly discuss the model (OQ) in the limiting case of strong rung 
coupling J± ^> J. In this limit the model (1) can be mapped onto the effective spin- 
chain Hamiltonian [H [16] . This allows to outline the symmetry aspects of the problem 
under consideration. 

At J± 3> J it is convenient to discuss the model by representing the site-spin 
algebra in terms of on-bond-spin operators [32] • Indeed an individual rung may be in a 
singlet or a triplet state with corresponding spectrum given by 



At h ext < J±, one component of the triplet becomes closer to the singlet ground state 
such that for a strong enough magnetic field we have a situation that the singlet and 



FigOb). 
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Figure 2. Qualitative sketch of the weak-coupling phase diagram of a two- leg 
ferromagnetic ladder as a function of the applied magnetic field (h ext ) and the isotropic 
interleg coupling (Jj_) [29 . In the case of isotropic spin-exchange interaction in legs 
(A = 1) the spin gapped Rung-Singlet (RS) phase is separated from the fully polarized 
ferromagnetic (FM) phase via the gapless Spin-Luttinger-Liquid (SLL) phase. In the 
case of finite intraleg exchange anisotropy (A / 1) the RS phase is splitted from the 
FM phase by the gapped Stripe-Ferromagnetic (SFM) phase. 



the S z = — 1 component of triplet create a new effective spin r = 1/2 system. One can 
easily project the original ladder Hamiltonian (pQ) on the new singlet-triplet subspace 

lt) = |s) = ^[|U)-|IT)], |^)s|r> = |u>- 

This leads to the definition of the effective spin 1/2 operators 

S+^ifl = (-l) n+a ^r+ , S^ a=h2 = i[7 + 2<] (4) 

The effective Hamiltonian in terms of the effective spin operators up to the accuracy 
of an irrelevant constant becomes the Hamiltonian of the spin 1/2 fully anisotropic XY Z 
chain in an effective magnetic field 

H eff= -^D^ r «+i + + Ar n' r «+i] + he// E r - ( 5 ) 

i n 

where h e ^ = h ext — J± + J/2. Note that in deriving we have used the rotation in 
the effective spin space which interchanges the x and z axes. 

At A = 1, the effective problem reduces to the theory of the XX Z chain with 
a fixed ferromagnetic XY anisotropy of 1/2 in a magnetic field, which allows for 
rigorous analysis [29]. The gapped phase at h e ^ < h e c {^ for the ladder corresponds 
to the negatively saturated magnetization phase for the effective spin chain, whereas 
the massless phase for the ladder corresponds to the finite magnetization phase of the 
effective spin-1/2 chain. The critical field h^ corresponds to the fully magnetized phase 
of the effective spin chain where the ladder is totally magnetized. From the exact ground 
state phase diagram of the anisotropic XX Z chain in a magnetic field |43j . it is easy 
to check that h^{{ 2 = =fJ/2. For the isotropic ferromagnetic ladder in a magnetic field 
this corresponds to a transitions from rung-dimers to SLL phase at h*f* = J± — J and a 
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transition from SLL phase into the fully polarized phase at h^f = J± (see FigfS]a). In 
this paper we compare this values of critical fields with the numerical results of exact 
diagonalization method. 

Away from the isotropic point A = 1 the effective Hamiltonian (jSj) describes 
the fully anisotropic ferromagnetic XYZ, chain in a magnetic field that is directed 
perpendicular to the easy axis. For the particular value of magnetic field h e ff = 0, the 
effective XYZ chain is long range ordered in Y — direction [H], accordingly the original 
ladder system being ordered in the direction perpendicular to the applied magnetic field 
with opposite magnetization on legs (stripe-ferromagnetic phase). For larger values of 
the effective field it is clear that this SFM order will be replaced either by the rung 
singlet phase or the phase with only one order parameter - magnetization along the 
applied field. 

The exact diagonalization study of finite ladder systems (up to length L = 12 (N = 
24)) show that the magnetic phase diagram of the ferromagnetic ladder with anisotropic 
legs in a transverse magnetic field really presents two quantum phase transitions with 
increasing magnetic field. The first transition is from a gapped RS phase to the gapped 
SFM phase and the second one from a SFM phase into the fully polarized FM phase 
(see Figfjb). 

3. Exact diagonalization results 

In order to explore the nature of the spectrum and the phase transition, we use the 
modified Lanczos method [4T1 14"5] to diagonalize numerically finite (N = 12, 16, 20, 24 
sites) ladder systems. The energies of the few lowest eigenstates were obtained for chains 
with periodic boundary conditions. 




Figure 3. Difference between the energy of the two lowest levels and the ground state 
energy as a function of the magnetic field, (a) h ext for J± = 2 J and A = 1.0 and (b) 
for J±_ = 3 J and A = 0.5 including different ladder lengths N = 12, 16, 20. 



To show, that the transition lines can be easily observed from the numerical 
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calculations of small systems we start our consideration by the case of isotropic chains. 
In this case, the magnetic field term commutes with the spin-spin exchange terms which 
makes the total x-component (Sf ot ) of spin a conserved quantity. Thus, the energy 
levels of ladder are labeled by S x ot . We have computed the vlaue of energy levels in the 
absence of magnetic field and added the effect of magnetic field by — Jh ext Sf ot . This 
makes very high accuracy for the value of energy levels. In Figj3^ we have plotted the 
three lowest levels of N = 12, 16, 20 ladder with J_l = 2 J as a function of the external 
field h ext . We determine the excitation gap in the system as the difference between 
the first excited state and the ground state. As it is clearly seen from this figure in 
the case of zero magnetic field the spectrum of the model is gapped. For h ext ^ the 
gap decreases linearly with h ext and vanishes at the critical field, h^f*. This is the first 
level crossing between the ground state energy and the first excited state one. To get 
an accurate estimate of h^f* we have obtained the first level crossing for system sizes of 
N = 12, 16, 20, 24. The finite size behaviour of this values leads us to h^* = 1.2±0.01 for 
N — > oo. The spectrum remains gapless for h^f* < h ext < h^f* and becomes once again 
gapped for h ext > h^f* = 2.0. We got h^f* = 2.0 as an exact value since there was no finite 
size correction at this value. We should mention that the critical field values obtained 
in the previous section come from the first order (perturbation) effective Hamiltonian 
approach and are not exact values. However, the values of the critical fields h^* and 
hgf* obtained from studies of the finite chains are very close to the values get in the 
previous section. With increasing field, for h ext ^> J±, J the gap increases linearly with 

It can be seen in Figj3^i for the gapless sector of the phase diagram at \f c xt < 
h ext < hgf* the two lowest states cross each other N/2 times, that the first crossing 
occurs at point h^f* and the last crossing occurs at the point h^f*. In this region 
we also observe numerous additional level crossing between the lowest second and 
third eigenstates. These level crossings lead to incommensurate effects that manifest 
themselves in the oscillatory behavior of the spin correlation functions. All crossings 
disappear at h ext > h. e c xt and the correlation functions do not contain oscillatory terms 
in this region of the phase diagram. 

In marked contrast with the isotropic case, the similar analysis of the few lowest 
levels for an anisotropic ladder in the presence of transverse magnetic field reveal a 
principally different behaviour. The energy difference between the two lowest levels and 
the ground state energy as a function of the magnetic field h ext has been computed for 
J_l = 3 J and J = 1.0 for different ladder lengths N = 12, 16, 20 and different values of 
the anisotropy parameter A = 0.3, 0.5, 0.7. As an example, in FigJSb we have plotted 
results of this calculations for A = 0.5. As it is seen from the figure, the excitation 
spectrum in this case is gapfull except at the two critical points h^* = 2.3 ±0.1 and 
n c2* = 3.0 ± 0.1. In the region of magnetic fields h^f* < h ext < h^f* the two lowest 
states form a twofold degenerate ground state in the thermodynamic limit and the 
spin gap, which appear at h ext > h^f*, first increase vs external field and after passing a 
maximum decreases to vanishe at h^f*. At h ext > h^f* the gap once again opens and, for a 
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sufficiently large field becomes proportional to h ext . These results are in good agreement 
with the results obtained in the studies of the fully anisotropic antiferromagnetic XYZ 
chain in a magnetic field |40j . 

To determine the properties of this model in different sectors of the phase diagram 
we have implemented the modified Lanczos algorithm of finite size ladders (N = 
12,16,20,24) to calculate the magnetizations M x,y,z and spin structure factors S xx (q), 
S yy (q) and S zz (q). 




Figure 4. a. The transverse magnetization M x as a function of applied field h ext for 
L = 10 ladder with J± = 3.0J and for different values of the anisotropy parameter 
A = 0.3, 0.5, 0.7. b. The intralcg spin structure factor in the x direction at q = for 
J± = 3.0J and A = 0.5 plotted as a function of the inverse ladder length L ( L = N/2 
) for different strength of the applied field h ext . The dashed line which corresponds to 
the h ext = h e c xt = 2.3J marks the transition from the RS phase to the intermediate 
phase with finite magnetization along the field and the dashed-dottcd line corresponds 
to h ext = h e c xt — 3.0 J and marks the transition into the phase with full magnetization 
along the field. 



In FigHh we have plotted the magnetization along the applied transverse magnetic 
field, M x vs h ext for N = 20, J± = 3 J and for different values of the anisotropy parameter 
A = 0.3,0.5,0.7. Due to the profound effect of quantum fluctuations the transverse 
magnetization remains small but finite for < h ext < \f c xt and reaches zero at h ext = 0. 
This is in agreement with results obtained within the weak-coupling continuum-limit 
bosonization treatment [29] and with the results for the magnetization obtained in the 
case of fully anisotropic XYZ chain [40J. For h ext > h^f* the magnetization increases 
linearly with increasing field once again in complete agreement with the predictions of 
bosonization analysis [29]. This behavior is in agreement with expectations, based on 
the general statement that in the gapped rung-singlet phase, magnetization along the 
applied field appears only at a finite critical value of the magnetic field equal to the spin 
gap [13-16]. However, in finite systems we do not observe a sharp transition at this point 
(h^f*) or close to the saturation value which happens at h ext > \f c xt . Magnetization along 
the directions perpendicular to the applied field remains small, but finite. However, as 
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Figure 5. The intraleg spin structure factor Sf^ g (q = 0) plotted as a function of the 
applied field h ext for L = 10 ladder with J± = 3.0J and for different values of the 
anisotropy parametere A = 0.3, 0.5, 0.7. 



we will show below, for Yi e c xt < h ext < h^f* the intraleg magnetization along the "y" 
direction is finite and shows the long-range ferromagnetic order. 

An additional insight into the nature of different phases can be obtained by studying 
the spin-spin correlation functions. In particular we study the magnetic field dependence 
of the different spin structure factors. We have calculated the intraleg spin structure 
factors defined by 



The field dependence of the intraleg spin structure factor S xx g (q = 0) is qualitatively 
the same as of the transverse magnetization M x . In Figfjb we have plotted the intraleg 
spin structure factor S xx g (q = 0) for different strength of the applied magnetic field 
chosen in the vicinity of h^f * c2 as a function of the ladder length L. As it is seen from 
this figure the dashed line which corresponds to the h ext = h^f 1 = 2.3 ±0.1 marks the 
transition from the RS phase to the intermediate phase with finite magnetization along 
the field and the dash-dotted line which corresponds to h ext = h e c xt = 3.0 ±0.1 marks 
the transition into the phase with saturate magnetization along the field. 

In Figj5] we have plotted S^ g (q = 0) as a function of the applied field h ext for 
L = 10, J± = 3J and different values of the anisotropy parameter A = 0.3,0.5,0.7. 
As it is clearly seen from this figure, in complete agreement with the bosonization 
results [29], there is not any long range ferromagnetic order along the "y" direction at 
h ext < h^f* = 2.3 ± 0.1 and h ext > h^' = 3.0 ± 0.1. However, in the intermediate region 
hpf* < h ext < h^2*, the spins of each leg show a profound ferromagnetic order in the "y" 
direction. 

To display the magnetic phase diagram of the anisotropic ferromagnetic ladder 




(6) 



ii, r 
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in a transverse magnetic field we have also calculated numerically the intrarung spin 
correlation functions. We have computed the on-rung dimerization order parameter 
given by 

dr = ^^2{0\ An" ( 7 ) 

n 

In Fig. [6^ we have plotted d r as a function of h ext for the ladder with J± = 3 J, A = 0.5 
and for different values of the ladder lengths N = 12, 16, 20, 24. As it is clearly seen 
from this figure for h ext < h e c xt , d r is close to —0.75) and the ladder is in the rung 
singlet phase. For h ext > h^, d r is slightly less then the saturation value d r ~ 0.25 and 
the ferromagnetic long range order along the "x" axis is present. Deviation from the 
saturation values —0.75 and 0.25 is the result of quantum fluctuations. In this case of 
finite systems and with chosen values of the rung exchange, the critical magnetic fields 
are high which strongly suppresses the quantum fluctuations. As a result the obtained 
averages of spin correlation functions are very close to their classical saturation values. 
In the intermediate range of field the system smoothly evaluates from a RS phase into 
the FM phase via formation of the stripe-ferromagnetic order. 




Figure 6. a. The rung-singlct order parameter d r as a function of the applied magnetic 
field h ext for the ladder with J± = 3.0 J, A = 0.5 and for different values of the ladder 
length N = 12,16,20,24. b. The rung-singlct order parameter d x , d\ and d z r as a 
function of the applied magnetic field h ext for the ladder with J± = 3.0 J, A = 0.5 and 
N = 24. 



To get additional data about the character of the spin ordering in the intermediate 
SFM phase we have calculated the "x","y" and " z" components of the rung-singlet 
order parameter given by 

d r = °l S l.n4n 1° > > J = *,y,Z. (8) 

n 

In FiglBb we have plotted d x , d y and d z as a function of the applied field h ext for 
N = 24 ladder with Jj_ = 3 J and for A = 0.5. As it is seen from this figure in 
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the RS, at h ext < h^f* the system clearly shows the SU(2) invariant order within the 
rung d x = d y = d z ~ —0.25. For h ext > h^ xt the SU(2) invariance is broken, d x 
increases with increasing field faster then the other two components and soon approaches 
the positive values corresponding to the ferromagnetic order in the "x" direction, 
while d y = d z remain negative indicating the antiferromagnetic correlations in these 
components. This shows the suppressed ferromagnetic order in the " z" direction and 
profound ferromagnetic response of the "y" direction presents completely convincing 
arguments in support of statement, that for intermediate values of applied field, at 
^ext < foext < ^ext f erroma g ne tic ordering in the "x" direction of spin from both legs is 
accompanied by antiferromagnetic order in the "y" direction of spins located on different 
legs and completely suppressed correlations in the " z" direction. It corresponds to 
the stripe ferromagnetic phase. Finally for h ext > h^f* the ordering in "y" and " z" 
directions are completely suppressed and the system shows ferromagnetic order with 
nominal magnetization per spin along the applied field. 

4. Conclusions 

We have studied the ground state phase diagram of the two-leg spin ladders with 
ferromagnetic anisotropic legs in a transverse magnetic field. We have implemented 
the modified Lanczos method to get the excited state energies as the same accuracy of 
the ground state one. Two quantum phase transitions in the ground state of the system 
with increasing magnetic field have been identified. The first transition is the gapped 
rung-singlet to the gapped stripe-ferromagnetic phase. The second - to the transition 
from the gapped stripe-ferromagnetic phase into the fully polarized ferromagnetic phase. 
This results are in complete agreement with the results obtained from the bosonization 
treatment [29] . 
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